Fluctuation persistent current in small superconducting rings 
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We extend previous theoretical studies of the contribution of fluctuating Cooper pairs to the 
persistent current in superconducting rings subjected to a magnetic field. For sufficiently small rings, 
in which the coherence length ^ exceeds the radius R, mean fleld theory predicts the emergence of a 
flux-tuned quantum critical point separating metallic and superconducting phases near half-integer 
flux through the ring. For larger rings with R^ the transition temperature is periodically reduced, 
but superconductivity prevails at very low temperatures. We calculate the fluctuation persistent 
current in different regions of the metallic phase for both types of rings. Particular attention is 
devoted to the interplay of the angular momentum modes of the fluctuating order parameter fleld. 
We discuss the possibility of using a combination of different pair-breaking mechanisms to simplify 
the observation of the flux-tuned transition in rings with > R. 
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I. INTRODUCTION 

The study of superconducting fluctuations has already 
a long history, for a comprehensive review see Ref. llj 
When approaching the superconducting phase from the 
metallic side, for example by lowering the temperature 
T, precursors of superconductivity reveal themselves long 
before the superconducting state is fully established. In 
this regime, electrons form Cooper pairs only for a lim- 
ited time. Being charged objects themselves, the Cooper 
pairs participate in charge transport. At the same time 
the density of states of the unpaired electrons is reduced. 
These simple qualitative arguments already indicate that 
superconducting fluctuations can affect both transport 
and thermodynamic properties of the metal outside the 
superconducting phase. Detailed studies of these effects 
have been conducted in different contexts. It is well 
known, for example, that fluctuation effects are more 
pronounced when the effective dimensionality of the su- 
perconductor is reduced or in the presence of disorder. 

In bulk superconductors the transition temperature Tc 
can be partially or even completely suppressed by var- 
ious pair-breaking mechanisms, most notably by apply- 
ing a magnetic field or introducing magnetic impurities. 
An additional pair-breaking mechanism can become ef- 
fective in doubly connected superconductors like super- 
conducting rings or cylinders, when they are threaded by 
a magnetic flux (p. In this case one observes so-called 
Little-Parks oscillations,'^ the transition temperature Tc 
is periodically reduced as a function of 0. Due to the 
periodicity, it is immediately evident that this effect is 
qualitatively different from the mere suppression of su- 
perconductivity by a magnetic field in bulk supercon- 
ductors. The period of the oscillations is equal to 1 as a 
function of the reduced fiux tp — (j) / (j)Q , where the super- 
conducting flux quantum is (po — n/e^ see Fig. [l] The 
maximal Tc reduction occurs when if takes half-integer 
values. 



The magnitude of the Tc reduction is size-dependent. 
It is convenient to measure the ring radius R in units 
of the zero-temperature coherence length ^ and to define 
r = R/^. A representative mean field phase diagram is 
displayed in Fig. ^ for two rings of different size. As 
we see in Fig. [T] mean field theory predicts a moder- 
ate Tc reduction for moderately small rings with r > 1. 
Most strikingly, it also shows that for very small rings 
or cylinders with r < 0.6 the transition temperature is 
expected to be equal to zero in a finite interval close to 
half-integer fiuxes [This regime is sometimes called the 
destructive regime.] Correspondingly, a flux-tuned quan- 
tum phase transition is expected to occur in these rings 
or cylinders at a critical flux ipco- The mean fleld tran- 
sition line can be found from Eq. (3.4) to be discussed 
below. 

As is well known, superconducting rings threaded by 
a magnetic flux (/) support a dissipationless persistent 
current.^ ^ In this article, we study theoretically the fluc- 
tuation persistent current in different regions of the phase 
diagram, both for rings with r > 1 and moderate Tc re- 
duction and for rings with r < 0.6 with strong Tc suppres- 
sion. In particular, we will study in detail the large fluc- 
tuation persistent current / which occurs even at fluxes 
for which Tc is reduced to zero while the system has a 
finite resistance. A short account of the most important 
results of this study has already been presented in Ref. |6l 
Here we extend our study to different regions of the phase 
diagram, discuss the results in a broader context, and in- 
clude details of the derivations. 

This study is motivated by recent experiments that 
are significant to our understanding of fiuctuation phe- 
nomena in superconductors with doubly-connected ge- 
ometry. Koshnick et al.'^ measured the persistent cur- 
rent in small superconducting rings with r > 1, for the 
smallest rings under study Tc was reduced by approxi- 
mately 6%. Superconducting fluctuations in rings of this 
size are by now well understood, both experimentally 
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FIG. 1: Mean field phase diagram. T^ip separates the metaUic 
(high r) and the superconducting (low T) phase as a function 
of the flux ifi = 4>/4>o through the ring. The transition line 



is determined by the condition Cqq = 0, cf. Eq. (3.2 1. The 



superconducting phase for small rings with effective radius 
r = R/^ < 0.6 is shown in dark gray. Mean field theory pre- 
dicts a full reduction of Tc for fiuxes between ipco ~ 0.83r and 
1 — ifico near (f = 1/2. In this article we focus on fluctuations 
in the normal phase for these small rings. The superconduct- 
ing phase for larger rings with r > 1 is shown in light gray. Tc 
is periodically reduced as a function of the flux ip, but super- 
conductivity prevails at low temp eratures. The dotted lines 
gives defined below Eq. for n G {0, 1, 2}. For r > 1 



it is well approximated by the formula ~ 7co(l — 'P /r ). 
The phase diagram is periodic in ip with period 1 for vanishing 
ring thickness. 



and theoretically.'^ll^ This is not so for smaller rings with 
r < 0.6. Strong Little-Parks oscillations for cylinders 
with r < 0.6, for which Tc is reduced to zero near half- 
integer flux, have been observed in a transport measure- 
ment on superconducting cylinders.'^^ It has so far, how- 
ever, not been possible to measure the persistent current 
in superconducting rings close to the flux-tuned quantum 
critical point. 

The difficulty to access the destructive regime for rings 
is that the experiments require a high sensitivity of the 
measurement device as well as low temperatures. At the 
same time, the magnetic field should be strong enough 
to produce a sufficiently large flux penetrating the small 
ring area. We address this issue in this manuscript by 
discussing the possibility that a combination of different 
pair-breaking mechanisms can lead to progress in this di- 
rection. More specifically, we consider the combined ef- 
fects of the magnetic fiux through the ring's center on the 
one hand and of magnetic impurities and/or the magnetic 
field passing through the bulk material of the ring on the 
other hand. By direct calculation, we further explore how 
the presence of the quantum critical point influences the 
persistent current away from the quantum critical point, 
e.g., for temperatures of the order of T° = Tc{(p — 0). 

The literature on superconductivity in systems with 
doubly connected geometry is extensive. We would like 
to point out a number of works, where related phenom- 
ena have been discussed. The possibility of finding com- 



plete suppression of superconductivity near half integer 
fiux in sm all su perconducting rings was pointed out by 
de Gennes.fii^ The phase diagram of superconducting 
cylinders was considered in Ref. taking into account 
the interplay of pair-breaking effects caused by the fiux 
on the one hand and the magnetic field penetrating the 
walls (of finite width) on the other hand. We will use 
their results when we discuss the influence of finite- width 
effects on the phase diagram. 

A detailed study of the fiuctuation persistent current in 
rings for which Tc is reduced to zero by magnetic impuri- 
ties (at any value of the flux) can be found in Refs.lMland 
1151 These works address a long standing puzzle related 
to the observation of an unexpectedly large persistent 
current in copper rings.'^It is suggested that these rings 
contain a finite amount of magnetic impurities, which 
suppress superconductivity and cause the rings to remain 
in the normal state even at low temperatures. Denoting 
the scattering rate on the magnetic impurities by I/tj, 
there is a critical rate I/t^c and an associated quantum 
critical point that separates the superconducting from 
the normal phase. If the measurements are performed 
on rings with a scattering rate that is larger but close 
to I/t^c; the corresponding fluctuations can lead to large 
currents in the rings. In contrast, in the parallel work 
of Ref. El and in the present manuscript we consider the 
opposite case, in which the phase transition is primarily 
tuned by the magnetic flux, so that for vanishing fiux 
and low temperatures the ring is in the superconducting 
state. We examine the influence of additional weak pair- 
breaking effects on the phase diagram and how they can 
help to experimentally observe the fiux-tuned quantum 
phase transition. 

In the experiment of Ref. 10 on cylinders it was ob- 
served that near half-integer fiux the resistance R along 
the cylinder drops as T decreases and then saturates for 
the lowest temperatures. In a later experimenl)^, regu- 
lar step-like features where additionally observed in the 
R — T diagram and interpreted as being due to a separa- 
tion into normal and superconducting regions along the 
cylinder. A number of theoretical works addresses the 
issue of transport in small superconducting-cylinders. In 
Refs. HSI and [THI the perturbative fluctuation contribu- 
tion to the conductivity of long superconducting cylin- 
ders near a flux-tuned quantum critical point was dis- 
cussed as a particular example for transport near a pair- 
breaking transition. In a broad sense, the general ap- 
proach is similar to ours, but the considered system has 
a different dimensionality and the work discusses trans- 
port, while we study a thermodynamic property. It has 
been suggested in Ref. \TE\ that the observable regime in 
the experiment!^ is dominated by thermal fluctuations 
and that at even lower T an upturn of R could be ex- 
pected. To the best of our knowledge, so far no detailed 
comparison between theory and experiment is available. 
The observed saturation in the experiment^*' corresponds 
to a strong reduction of the normal resistance and as such 
lies outside the region of validity of the perturbative ap- 
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proach. The role of inhomogeneities along the cylinder 
axis has been further emphasized in Ref. [20l In Ref. |2T] 
a mean field model was proposed that takes into account 
inhomogeneities caused by a variation of parameters like 
the mean free path or the width along the cylinder axis 
and a specific profile was found that would quantitatively 
fit the experimental phase diagram both as far as the sat- 
uration and the step-like features are concerned. Fluc- 
tuation effects, on the other hand, where neglected. It 
seems likely that a complete description would have to 
include both inhomogeneities and fiuctuations, which is 
very demanding. In this respect, the situation with rings 
is more advantageous. Since the typical size of inhomo- 
geneities is larger or of the order of the coherence length, 
they are unlikely to play a role for the superconduct- 
ing rings under study here and we may focus on fluc- 
tuation effects only. We make detailed predictions for 
flux and temperature dependence of the resulting fluc- 
tuation persistent current. As mentioned before, due to 
experimental difficulties, (to the best of our knowledge) 
no measurements of the destructive regime are available 
for rings yet. 

Fluctuation effects in moderately small superconduct- 
ing rings (r > 1) were studied in Refs. [9) and [22] with 
particular emphasis on the regime of strong fluctuations 
near the transition. These works introduce the idea that 
the strong fluctuation regime near the thermal transition 
can be described in terms of one or two coupled angu- 
lar momentum modes of the order parameter field in the 
classical GL functional. Wc will make use of this idea 
and derive more detailed results for the persistent cur- 
rent and susceptibility close to integer and half-integer 
fiuxes. 

Ref. [S] studies superconducting fiuctuations in rings 
near the thermal transition using a numerical approach 
based on the mapping of the classical GL-functional 
onto the problem of solving an effective Schrodinger 
equatiorp2l. The results of this approach agree well with 
the experiment of Koshnick et al™ This approach works 
nicely in cases where many angular momentum modes 
of the order parameter field give a sizeable contribution 
to the persistent current, but can become cumbersome 
in the opposite limit (see Ref. In this sense, it is 

complementary to the approach used in this manuscript 
(as well as in Refs. I6l9l22p . which is well suited for rings 
that are so small that only one or two angular momentum 
modes are important. 

The role of the back-action effects caused by the self- 
induction of cylinders and rings was discussed in this 
context in Refs. [55] and 26. While we will discuss finite 
thickness effects, we will generally assume that the self- 
induction and the persistent currents of the rings under 
study in this article are sufficiently small so that back- 
action effects can be neglected. 

The article is organized as follows. Section |ll] is de- 
voted to the persistent current in rings with r > 1. Some 
calculational details are presented in appendix |X] In sec- 
tion III we study in detail the fluctuation persistent cur- 



rent in different regions of the phase diagram for rings 
with r < 0.6, including the vicinity of the quantum crit- 
ical point. Details of the calculation are relegated to 
appendix |B] In section |IV| we discuss different ways to 
reduce Tc to zero in rings with finite width or by intro- 
ducing magnetic impurities. 



II. THERMAL TRANSITION FOR RINGS 
WITH r > 1 

In this section we will discuss the description of rings 
with only a moderate suppression of Tc, i-e., rings for 
which r > 1. A condensed discussion has already been 
presented in Ref. 'W. Here we take the opportunity to 
provide additional information. For the rings with r > 1 
the superconducting transition occurs at a finite temper- 
ature and fiuctuations can be described with the help of 
the classical GL functional in which the order parame- 
ter field is static. In the imaginary time formalism this 
amounts to neglecting order parameter field components 
with finite Matsubara frequency in the functional- (a 
formal justification will be given in the first paragraph 
of Sec. Ill C| below). This simplified description is valid 
close to the transition line in the Tc — (p phase diagram. 



A. Ginzburg-Landau functional 

The starting point for our discussion of superconduct- 
ing fluctuations in rings with r > 1 is the classical GL 
functionaP^l 



Z = j D{^,'tP*)exp{~T/T) 



(2.1) 



T ^ Tn + dr (a\ij{r)\ 



1 

Am 



(-iV-2eA(r))V'(r)| 



describes the the normal (non-superconducting) part 
of the free energy. It gives rise to a normal component of 
the persistent current.^* Since we are mainly interested 
in a regime close to the transition, however, the fluctua- 
tion contribution is much larger and we will not discuss 
J-jy further. At the mean field level, the sign change of 
the quadratic form in ip signals the onset of the super- 
conducting phase, which motivates the parametrization 



aT^s, where e 



A characteristic length 



scale, the (zero temperature) coherence length, can be 
identified £_ = l/y^AmaT°. The microscopic theory for 
disordered superconductorJ^SI gives rise to the following 
relations 



(2.2) 



where D is the diffusion coefficient. The normalization 
of Ip allows for a certain arbitrariness, this is why only 
the ratio of and b is fixed. 
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The quartic part of the functional stabihzes the system 
once it is tuned below the transition temperature. Above 
this temperature, the quartic term gives only a small con- 
tribution to thermal averages, except in the very vicinity 
of the transition, the so-called Ginzburg region. As long 
as one stays outside of this region on the metallic side, 
one can restrict oneself to a quadratic (Gaussian) theory 
(i.e. neglect the quartic term), which is much easier to 
handle theoretically, but becomes unreliably close to the 
transition where the quadratic theory becomes unstable. 
Importantly, even above the transition temperature, in 
the normal region of the mean field phase diagram, the 
average of |V'P with respect the functional T is finite. 

After this preparation, we turn to the description of 
superconducting rings. When the superconducting co- 
herence length ^(T) — (,l\fe and the magnetic pene- 
tration depth A(r) are much larger than the ring thick- 
ness, the system is well described by a one-dimensional 
order parameter field ';/': albeit with periodic boundary 
conditionJ^. In order to account for these boundary con- 
ditions, it is convenient to introduce angular momentum 
modes as 



(2.3) 



V — 2ttRS±^ is the volume of the ring, the cross- 
section of the wire forming the ring. The vector potential 
can be chosen as A = B x r/2, the integration as / dr — > 
S'j^i? Then, the free energy functional takes the 

form 



2V ^ 

nmkl 



-k.l+mlpn^l^^Tpk^l^! (2.4) 



where a„^ = a + (n — ip)'^ /2mR^ . Let us make three 
important observations. First, the free energy functional 
is flux dependent. As a consequence a persistent current 
can flow in the ring. Second, the functional is periodic in 
the reduced flux (p with period 1. Correspondingly, the 
same is true for all thermodynamic quantities derived 
from the the GL functional. This property holds strictly 
speaking only in the idealized limit of a one-dimensional 
ring. In reality, the external magnetic field also pen- 
etrates the superconductor and provides an additional 
mechanism for the suppression of superconductivity. We 
will come back later to this point in section [rv] The third 
observation is that now the kinetic energy of the Cooper 
pairs vanishes only when the reduced flux ip takes inte- 
ger values. Otherwise it gives a finite contribution to the 
quadratic part of the functional and therefore the transi- 
tion takes place at a temperature Tc{ip) that is in general 
reduced with respect to of the bulk material. Let us 
parameterize 

and a„ change sign at a temperature 



aT^En^- Then = {T - Tn^)/T^ 



r°[l-(n-^)Vr^ 



(2.5) 



This temperature can loosely be interpreted as the tran- 
sition temperature of the nth angular momentum mode 



■>pn- The mean field transition for the ring occurs at 
Tcip that is equal to the maximal r„ for given tp, i.e., at 
the point where the first mode becomes superconducting 
when lowering the temperature (cf. Fig. [T]). 

In the subsequent discussion, the parameter A = 
1/r^Gi will play a crucial role. Its relevance is now eas- 
ily understood. is a measure for the typical spac- 
ing between the transition temperatures T„ for differ- 
ent modes, since e.g. [To{(p) - Ti{p)]/T^ = (1 - 2p)/r'^ 
(compare Fig. [I]). This spacing can be compared to the 
typical width of the non-Gaussian fluctuation region, Gi. 
Only in this region fluctuations are strong. The zero- 
dimensional Ginzburg parameter Gi of relevance here is 
given 



Gi 



2b 



^ 7C(3) 



(2.6) 



where v is the density of states at the Fermi level and 
V is the volume of the ring. The parameter A = 1/r^Gi 
determines whether a description in terms of a few modes 
only is a good approximation in the critical regime (for 
A ^ 1) or not. Defining the dimensionless conductance 
of the ring as 5 = Rq/Ro = 2e^ v D RqV / {27r Rf , Rq = 
Tr/e^, one can find an alternative expression, A w 5y^/r. 
We state this alternative (but equivalent) expression for 
A, because it might be more convenient for estimates 
when performing an experiment. 

An analytic computation of the functional integral nec- 
essary to obtain Z is in general not possible and one has 
to resort to approximation schemes. If the spacing is 
large (A ^ 1) and one is interested in the region of strong 
fiuctuations close to the transition temperature Tc{p), an 
effective theory including only one angular momentum 
mode ■(/'„ for n ^ p is applicable. 



^ n 



2V 



(2.7) 



This is so, since in this case the temperature T ~ Tc{p) 
lies far above the individual transition temperatures 
Tm{'^) (jn ^ n) of all other modes ipm and they will give 
only a small contribution when calculating observables. 
This is very convenient, because in this case one comes 
to the zero-dimensional limit of the GL functional .1^^ 
The partition function based on this free energy func- 
tional can be calculated exactly and all thermodynamic 
quantities derived from it. 

Clearly, for half-integer values of the flux the spacing 
between two adjacent modes always goes to zero and due 
to this degeneracy at least two modes are required for the 
description. Let us choose for definiteness the example 
< p < 1, then one may work with 



-^01 



1=0,1 



In this situation the parameter w Xjhr^Gi is still 

useful, because if it is large, additional modes need not be 
taken into account and a two-mode description is valid. 
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If ^/g/r is not exceedingly large, the small contribution 
of the remaining modes can easily be accounted for by us- 
ing the Gaussian approximation for them. One should, 
however, not forget that the presence of the dominant 
mode(s) can influence the effective transition tempera- 
tures of all others via the quartic term. As an example, 
let us write the resulting effective action for the modes 
with n 7^ assuming that ~ and the mode with 
n = is the dominant one. 



(2.9) 



Essentially the same argument was first presented in 
Ref. m If temperatures are sufficiently high T >• T^{\ + 
Gi) the quartic term may be dropped altogether and 
one may work with a purely Gaussian theory Tn ~ 



B. Persistent current and susceptibility 

The persistent current / is found from the free energy 
F — —T\nZ by differentiation / = —dF/dcf). The nor- 
malized current is given by 



I = I/{T^J<t>,) 



2a 



(2.10) 



The averaging is performed with respect to the functional 
F in Eq. (2.4). Just as the free energy functional F^ 
the persistent current i is periodic in the flux with 
period one. Since it is also an odd function of (/?, the 
persistent current vanishes when ip takes integer or half- 
integer values. 

a. Case Ki n, T ^ Tc'. As pointed out above, 
the most important contribution in the regime of non- 
Gaussian fluctuations close to integer fluxes comes from 
the angular momentum mode il^n with the highest tran- 
sition temperature Tn^p. One may then approximate 



Eq. (2.4) by a single-mode and calculate with Fn = 
an\i'nr+ This is the Od limit of the GL 
functionapSl already introduced above. In this limit. 



Eq. (2.101 gives 



(2.11) 



(2.12) 



in = 4A(n - ip)f{xn) for p n. 

Here a;„ = £n/Gi and the function 

^ exp(-x^) ^ ^ 
-y7rerfc(a;) 

is defined with the help of the conjugated error 
function.'^ All persistent current measurements will fall 
on the same curve, if the persistent current - measured 
in suitable units i = I/{T^/4>q) - and the reduced tem- 
perature — {T — Tcip)/T^ are scaled as 

T 

i^i — , e^^e^ry^. (2.13) 

This relation can serve as a valuable guide in character- 
izing different rings in experiments. 



b. Gaussian theory for T ^ , estimate for T <^ 
r° ; It is possible to make contact with the Gaussian and 
the mean field results using the asymptotic expansion of 
the conjugated error function 

^/nxerfc{x) « exp{~x^){l - l/{2x^)) {x cx3X2.14) 

and the limit erfc(x) — >■ 2 for x — > — oo. Far above Tc 
one obtains as a limiting case the Gaussian result for 



a single mode 



2{n - ip)/r^ 



■nip 7 



that can also be 



obtained directly by neglecting the quartic term in the 
GL functional. In this form, however, it is of limited use, 
since for the temperatures in question one should sum the 
contribution of all modes. Indeed, in this case one can 
use the relation (iV'nP) ~ T^/o-n {o-n was defined below 
Eq. 2.4 1 and perform the sum in Eq. (| 2.10[ ) to obtain a 
result that is valid at arbitrary fluxeS^^i 



i{(p) 



-2tt sin(27r(/3) 



cosh(27r-yer) — cos(27r93) 



(2.15) 



One of the main features of this result besides the pe- 
riodicity in if is the exponential decay of the persistent 
current as a function of temperature for e > l/(27r)^r^, 
which is due to a mutual cancelation of the contributions 
of many modes to the persistent current. Turning back 
to Eq. (2.11), we see that far below Tc one recovers the 



mean field result 



r^Gi 



2 '^'^ 



(n-ip), 



(2.16) 



which gives an estimate for the persistent current in 
the superconducting regime. In this approximation the 
current grows linearly with \T — Tc\ and as soon as 
|e| Gi + {n — (p)'^/r'^ (i.e. |a;„| >• 1) one expects a 
sawtooth-like behavior as a function of the flux (i.e. lin- 
ear dependence from ip = n — 1/2 to (p = n + 1/2 passing 
through zero at integer n) with a discontinuous jump at 
half integer p. Both the Gaussian and the mean field 
result are reliable only outside the region of strong fiuc- 
tuations, the persistent current i„ in Eq. (2.11) covers 



this region and interpolates smoothly between them. 

c. Case ip K n + 1/2, T Tc: At half integer values 
of p, the transition temperatures for two modes become 
equal. In the vicinity of this point in the phase diagram 
the two dominant modes influence each other, their cou- 
pling becomes crucial. We discuss the case (p ^ 1/2 for 
definiteness, and use the form of the free energy func- 
tional already displayed in Eq. (2 



Calculation of the persistent current in the presence 
of the coupling requires a generalization of the approach 
used for the single mode caseW Explicit formulas for the 
persistent current are derived and displayed in appendix 
\K\ for the sake of completeness. For a graphical illustra- 
tion see Fig. 2 of Ref. "51 We define the susceptibility as 
X = —dl/dcf). Differentiating the expression for Z2 (see 
appendix [A|) one obtains 



Xtp=l/2 



X 



= 4A 31 (a;)~4A2 (a;), (2.17) 
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wher^ 

a; = a;o(l/2)=xi(l/2) = ;^ + ^A. (2.18) 

Note that x = at the mean field transition for ip = 1/2, 
i.e. the condition a; = defines Tc,i/2 = ^"(1 - l/Ar"^). 
The dimensionless smooth functions gn are defined by 
the relations 



51(2;) 
32(2;) 



1 



2J{x) 
3 



erfc(a;) 



2y/^J{x) 



2x 

" y 

3a; 
2J(x)* 



(2.19) 

erfc(a;) — 1, 

(2.20) 



where J{x) — dt es* erfc(t). Exact results can be 
given for the functions gn at the transition, i.e. for a; = 0. 
These give a useful estimate for the magnitude inside the 
fluctuation region, 51 (0) = -/7r/[2-\/3arctanh (l/\/3)] w 
0.78 and 32(0) = v^/[2arctanh (l/Vs)] - 1 ~ 0.315. For 
large A = 1/r^Gi « 5y^/r one can neglect the first term 



in Eq. (2.17). Then one obtains Xi/2 = ~ 4A^ ff2 (a;). 
For the susceptibility close to integer flux one easily ob- 
tains Xo = 4A/(a;o) from Eq. (2.11). Comparing to the 
expression for Xi/2: 



X1/2/X0 ~ -2.77g/r. 



(2.21) 



Experimentally, a strong enhancement of the magnetic 
susceptibility near y = 1 /2 compared to 1^9 ~ was 
observecPH and Eq. (2.21) demonstrates that it is con- 



trolled by the parameter ^Jgjr. If it is large, the current 
rapidly changes sign as a function of the flux at half- 
integer flux, leading to a saw-tooth like shape of i^p . The 



full T dependence of Xip=\ji is given in Eq. (2.17). The 
shape of the persistent current as a function ot the flux 
has been discussed in more detail in Ref. [6l Let us just 
stress the main physical mechanism at work here. To 
be specific, we discuss the vicinity of — 1/2. Close 
to 1 /2 there is a competition of two angular momentum 
modes, -00 and -01, that are almost degenerate. If one 
tunes to a slightly smaller fiux, say, then the mode -00 is 
dominant, because oq is smaller than a\. The effect of 
the coupling term in this case is to further weaken the 
mode V'l (this can be seen from an effective action in the 
form displayed in Eq. |2.9| ), for a repulsive interaction 
the dominant mode suppresses the subdominant mode. 
Since the contribution to the persistent current of these 
two modes is opposite in sign, the result is an almost 
saw-tooth like shape of the persistent current.t^. Similar 
effects related to the competition of two order parameter 
fields have been discussed in the past, see, e.g., Ref. [Ml 
d. Summary: To summarize, in this section we ar- 
gued that for A — hyfgjr ^ 1 one can use simple approx- 
imation schemes to calculate the persistent current and 
susceptibility near T^. The significance of the parame- 
ter A is as follows. If A is large, the statistical weight 
of one angular momentum mode in the strong fiuctua- 
tion region by far exceeds the weight of all other modes. 



unless two modes become degenerate. The degeneracy 
points are (j) = n + 1/2. Our calculations were based on 
formula [2 . 1 0| which is a general expression for the fluctu- 
ation persistent current expressed in terms of the ther- 
mal averages (j'i/'np). Near integer fiux {ip « n) and for 
T « Tc it is sufficient to include only one mode, which 
leads to formula |2.11| The periodicity of the phase dia- 
gram is not crucial here. Instead, the relevant free energy 
(Eq. 2.7) has the same form as that of a small supercon- 
ducting grain in the zero-dimensional limit .fSfil This is a 
drastic simplification compared to the original problem 
and implies a high degree of universality. With the ap- 
propriate scaling given in Eq. 2.13[ i((^)-data measured 



for rings with various parameters should fall on the same 
curve. As the temperature increases or the flux comes 
closer to the degeneracy points, the restriction to only 
one mode is no longer a good approximation. Formula 
|2.15[ which is valid for high temperatures, includes the 
Gaussian contribution of formula 12.111 and all the other 
modes. It is therefore applicable for arbitrary fluxes. The 
shape of i{ip) becomes more and more sinusodial as the 
temperature increases, this is the result of the combined 
contribution of many modes. Returning to the vicinity 
of Tc{(p), in formula 



2.17 



we calculate the contribution 
of two modes to the magnetic susceptibility at half in- 
teger fluxes {(fi n + 1/2). Formula 2.21 compares the 



magnetic susceptibility at half-integer and integer fluxes 
and reflects the large enhancement of the susceptibility 
near ip n + 1/2, which is controlled by two fluctuating 
angular momentum modes. We interpret this enhance- 
ment as being due to the suppression of the subdominant 
by the dominant mode as a result of their interaction. [M] 
This interaction is induced by the quartic term of the GL 
functional. In table |I] we summarize the expressions and 
analytical approximations for the fluctuation persistent 
current near the thermal transition together with their 
range of validity. 



III. GAUSSIAN FLUCTUATIONS FOR RINGS 
WITH r < 0.6 

So far we have discussed the limit of moderately small 
rings with r = R/S, > 1- For these rings the transition 
temperature is only weakly suppressed at flnite flux and 
the phase transition occurs at a finite temperature. We 
will now discuss even smaller rings with r < 1. For such 
rings the theoretical description based on the classical GL 
functional we used so far is not valid in large regions of 
the phase diagram, as it is applicable only in a relatively 
small temperature interval close to T". For rings with 
r < 1, however, the transition near half-integer fiux can 
occur at temperatures far below T^?, or even at vanishing 
temperatures. The theoretical description for these rings 
can be developed in close analogy to the general theory 
of pair-breaking transitions. In this section we will first 
determine the mean field transition line and then cal- 
culate the contribution of Gaussian fluctuations to the 
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Formula 


Meaning 


Range of validity 


(2.10) 


PC, general formula, all ij^n contribute 


thermal transition 


(2.11) 


PC, strong fluctuations, one ipn dominates 


ifxn, \T-T^\ <T,Gi,K-> 1 


(2.15) 


PC, gaussian fluctuations, all contribute 


tp and A arbitrary, T > (1 + Gj)Tc 


(2.16) 


PC, mean field estimate, one ■;/>„ dominates 


99 « n, r < T„ A > 1 


(A5) 


PC, strong fluctuations, two ip^ contribute 


ip^n + 1/2, |T - Tel < T^Gi, A » 1 


(2.17) 


Susceptibility, strong fluctuations, two contribute 


(p^n + 1/2, \T - r,| < T,Gi, A > 1 



TABLE I: An overview of different formulas for the fluctuation persistent current (PC) near the thermal transition for rings 
with r > 1. The range of validity for the analytical approximations is also estimated. Tc is the (flux dependent) critical 
temperature, ip = 4>/{h/2e) is the dimensionless flux through the interior of the ring and Gi is the Od-Ginzburg parameter 
defined in Eq. (2.6). The dimensionless parameter A has been introduced after Eq. (2.6) and can be much larger than 1 for 
small rings. 



persistent current outside the superconducting regime. 



C ^(0,0) ~ 0. It defines an implicit equation for Tc{tf) 



A. Mean field transition line 

The partition function is conveniently formulated in 
terms of an integral over a complex order parameter field 
A as Z = / D(A, A*)exp(-S'), where 

S = T^^A*(n,w)/:-i(n,w)A(n,a;). (3.1) 



The field A - unlike the field in the classical GL func- 
tional of Eq. |2.1| - is dynamical, = 27rmT is a bosonic 
Matsubara frequency. For a derivation of Eq. |3.1| see, e.g.. 
Chapter 6 of Rcf. 1, Let us note here that the field V'o 
is proportional to the static component A(w„ = 0)- A 
neglect of fields with nonzero Matsubara frequencies can 
be justified for the thermal transition, where it leads to 
the classical GL functional, but is not justified near the 
quantum phase transition. In full analogy to previous 
considerations, A has been expanded in terms of angular 
momentum modes. The fluctuation propagator C fulfills 



= ln^^ 



(3.2) 



1 



;|/2 



J ' " \2 ' 2tiT J ' \2, 

Here we introduced the pair-breaking parameter 

an = erin - ifif /2 (3.3) 

for the problem under consideration and v is the den- 
sity of states at the Fermi energy. We use the notation 
Et = D / and ?A is the Digamma function. Since we 
are interested in the mean field transition line and in the 
Gaussian fluctuations on the normal side of the transi- 
tio n, th e quartic and higher order terms in the action of 
Eq. 3.1 may be safely neglected (see, e.g., the booli^'). 

As usual, the mean field transition occurs when 
C^^{n,ui) changes sign first for arbitrary n and w. As- 
suming that |(^| < 0.5, this happens for n = 0, w = 0, 
so that the condition for the mean field transition reads 



(3.4) 



The larger ao, the lower values of Tc{(p) are required to 
fulfill this relation. Eventually, for 



ao(</'c) = 27rr,"exp(V(l/2)) = nT^j2^E 



(3.5) 



the transition temperature T^^if) vanishes, i.e. one 
reaches a (flux-tuned) quantum critical point. This hap- 
pens for the critical flujP 



Vc = 7rr/(2A 



7B « 1.78. 



(3.6) 



Due to the flux-periodicity of the phase diagram, a quan- 
tum transition can only be observed in the ring geometry 
if (pc < 1/2, which implies r < y^St^^/tt w 0.6 (compare 
Fig. [T]). Notice that this critical value of r = i?/^ is less 
restrictive than a naive application of the quadratic ap- 
proximation valid for r ^ 1 would suggest. The latter 
would give 1 - (l/2r)2 = Q =^ r = 1/2. 



B. Persistent Current 

The formula for the persistent current / 
in the Gaussian approximation reads 



dip 



(3.7) 



where we used Eq. |3.1| In order to gain a good qual- 
itative and quantitative understanding of the flux and 
temperature dependence, the magnitude and relevance of 
the scales involved in the problem, we will in the follow- 
ing derive simpler expressions for the persistent current 
and discuss various limiting cases. A particular empha- 
sis will be put on the analysis of the persistent current 
in the vicinity of the quantum critical point. This ap- 
proach will be corroborated by a direct numerical eval- 
uation of (3.7). When doing so, some care needs to be 
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exercised in order to correctly deal with the slow conver- 
gence properties for large values of \n\ and |a;|. As before, 
we find it convenient to consider the dimensionless quan- 
tity i — I/{T^/(Pq) and for definiteness analyze fluxes if 
in the interval (0,0.5). Since the persistent current is 
periodic i{(p + 1) = i{(p) and odd i{ip) = —i{~ip) in the 
flux, this is sufficient to infer the persistent current for 
arbitrary fluxes. 

As is shown in appendix |B] the persistent current can 
be written as the sum of two contributions, i = ig + ins- 
The rationale behind this decomposition is the follow- 
ing: The singular part is diverges on the transition line 
tpc{T). Outside the Ginzburg region on the normal side 
of the transition, is is still strongly flux and temperature 
dependent. It describes by far the dominant contribution 
to the fluctuation persistent current in the entire normal 
part of the phase diagram. The nonsingular part i„s, on 
the other hand, displays a smooth flux dependence and 
is much smaller, is and ins have opposite signs. We will 
mostly discuss is, for explicit formulas for ins we refer to 
appendix [B] The expression for is reads 



sin(27r <y9) 



cosh(27rz5) — cos(27r(^) 



(3.8) 



Here, Zg is defined as the unique positive solution of the 
equation 



In 



T 
TO 



AttT 



(3.9) 



for which the argument of the digamma function on the 
right hand side is positive. The formula for ig in Eq. |3.8| 
can be viewed as a generalization of Eq. |2.15| to the entire 
normal part of the phase diagram (outside the Ginzburg 
regime). 

In order to find a more convenient expression for Zg 
it is useful to define the function adT) describing the 
phase boundary in the aQ-T phase diagram, i.e. ac{T) — 
ao{(pc{T)), 



In 



T 



«c(T) 



27rr 



One immediately reads off 

Zs = Vd^l - 2«c(T))/eT, 



(3.10) 



(3.11) 



Here, adT) in Eq. ( 3.10 ) is allowed to become negative as 
soon as T > T°. Note that Zs is always real for T > T^. 
This is no longer true for T < T^. In this regime it is 
instructive to use the fact that adT) and the temper- 
ature dependent critical flux (fd"^) closely related, 
ETfliT) = 2adT). As a resuh, 

Zs = cpdT)V\Lo\/2adT)-l, {T<T^). (3.12) 

Now Zs is real for oj > 2adT), but becomes purely imagi- 
nary for uj < 2adT). The fact that Zs becomes imaginary 



for small |w| and T < T°, but not for larger T > T^, is 
intimately related to the occurrence of the phase transi- 
tion. Indeed, the denominator in the expression |3.8| for 
is vanishes for ip — (fd"^) at w = 0, signaling the onset 
of the superconducting regime. For reference, recall that 
a,o = adO) = ^T0/27b « 0.88TO. 



C. Fluctuation persistent current for T > 



Let us first make contact with the Gaussian result for 



larger rings r > 1 stated before in Eq. (2.15). For T 



when the logarithm on the left hand side of Eq. (3.10) 



7^0 



is small, one can expand the Digamma function on the 
right hand side. When keeping only the most dominant 
term in the sum, the term with vanishing Matsubara 
frequency, one easily finds Zs ^ ^fer and in this way 
reproduces the result of Eq. (2.15) after setting t ~ \. 
The restriction to a; 

Zs(uJn) ^ \f\y ~ 



= is justified in this case, because 
■Kr^/n/2 is real and larger than one 



for finite Matsubara frequencies, and the denominator in 
the expression for is becomes large, thereby strongly sup- 
pressing the contribution of finite n 7^ 0. 

Now we turn to the smaller rings with r < 0.6, for 
which a quantum phase transition takes place at zero 
temperature. We first analyze the flux dependence of i 
for a given temperature T ~> T^. As long as (pco ^^nd 
correspondingly r do not become very small, the same 
argument concerning the importance of the lo — Q com- 
ponent that was used for rings with r > 1 in the pre- 
vious paragraph is applicable here. For an estimate, 
Zs{loi) ~ ^J\loi\/£t = 7rr/2 becomes equal to 1/2 only for 
r < I/tt, which corresponds to </3c0 ~ 0.27. Interestingly, 
the same parameter ■Kr/2 determines the relevance of the 
nonsingular contribution in this case (see appendix [b| . 
As long as this parameter does not become considerably 
smaller than one, it is therefore safe to concentrate on 
the a; = term of the singular contribution only. When- 
ever it is justified to use only this term near T^, then it 
is also justified for larger temperatures (as can be seen 
by comparing Zs to zq defined in appendix [b|) . 

To summarize this somewhat technical discussion, for 
not too small rings with Lpco ^ 0.3 we obtain a very good 
description for the entire temperature range T > T^ -\- Gi 
by the formulcP2l 



-27ri- 



sin(27r(y9) 



h{2-K ^ -2adT) I £t) - cos(27r¥3) ' 



COS. 



(3.13) 



The formula given above remains valid for T < T^ for 
fluxes, for which the ring is in the normal regime. In 
Fig. ([2]) we display the flux dependence of the persistent 
current for different temperatures. At small fluxes, the 
persistent current is proportional to ip. For temperatures 
close to r°, this leads to a rapid increase of |i| for small 
Lp. As the flux increases, however, Td'p) decreases rapidly 
for the small rings under consideration here. Therefore 
the distance to the critical line in the phase diagram 
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FIG. 2: The modulus of t he fluctu ation persistent current 
i = I/[T°/cj>o] (see Eqs. 



3.7 



and 



3.21 is shown for a ring with 



ipco ~ 0.4, corresponding to r « 0.48, for three temperatures 
above T° [t = T/T°] (soHd lines). The da shed lines show the 
thermal {oj — 0) part of is, see Eq. ( |3.13 1, which gives a very 
good approximation. The current decreases as the tempera- 
ture grows. At the same time the shape changes considerably 
as explained in the text. In the inset the phase diagram is 
displayed, indicating the temperatures in question. 



grows as ip increases while the temperature is kept con- 
stant. This is why the current subsequently drops. For 
larger temperatures, the situation is different. As the 
temperature increases, cosh(27rZs) « exp(27rzs) grows 
and the flux dependence is determined by the numerator 
in the expression for is, Eq. ( |3.8[ ). As a consequence the 
shape becomes more sinusoidal. For an e stima te we can 
use that for T > r°, z,. - y/er, see Eq. ( |2.15[ ). The ex- 
ponential decay of i and transition to a sinusoidal shape 
therefore starts at e « l/(27rr)^. The persistent current 
at these high temperatures results from a combined ef- 
fect of many angular momentum modes and is therefore 
much less sensitive to the rapid drop of the transition line 
than for temperatures T ~ T^. 



D. Fluctuation persistent current near the 
quantum critical point and at intermediate 
temperatures for rings with r < 0.6 

The situation is quite different in the low temperature 
limit T <^ for rings with r < 0.6, to which we will turn 
now. Indeed, for very low temperatures a restriction to 
the thermal fluctuations, namely those with cj = 0, is not 
justified as we will see now. 

We discuss the vicinity of the critical line for T <C 
T°. Here one can expand the denominator in the general 
expression |3.8| for is in small 

AcpT^[p^iPciT)]/iPciT) (3.14) 

as well as small |a; | /etc (T). In this way one obtains the 
approximate relatiorP^I 



2pc{T) 
IepIq 



h{ApT,t), 



(3.15) 



where the dimensionless function h is defined as 

/i(Av3T,i) =7rry — ^ -r^—^^ . (3.16) 

The same expression can be obtained directly from the 
initial formula for i {i ^ I/{T^/(j)Q) with / given in 
Eq. 3.7), if one identifies p ~ (pc{T), and considers the 
most singular angular momentum mode n = only. It is 
worth noting, however, that for fixed n the sum in uj is ul- 
traviolet divergent (even before expanding in \Lo\/ac{T)). 
When proceeding in this way a cut-off has therefore to be 
introduced by hand. In contrast, our formula (Eq. 3.8) 
for is immediately reveals that terms in the sum with 
UJ > 2ac{T) are suppressed, since Zg becomes real and 
cosh(27r2;s) grows rapidly for larger Matsubara frequen- 
cies. We can therefore perform the sum with logarith- 
mic accuracy and choose to = 2ac(T) as the upper cut- 
off. Put in different words, the upper limit for the 
summation is effectively provided by mutual cancelations 
between different angular momentum modes. 
The result of the described procedure is 



h{AipT,t) 



1 

2^ 



+ tp[l + 



2A(pT 



where 



2Aipx OLc{T) AipT 



t 



ttTO 



V(l + s),(3.17) 



(3.18) 



The first term in the expression for h is the classical lu = 
contribution to the sum. These thermal fluctuations 
are proportional to the temperature and correspondingly 
vanish for T — > 0. This does, however, not imply the 
vanishing of the persistent current in this limit. In order 
to see this more clearly, let us display the asymptotic 
behavior of the function h: 



^ + In ^ 

1 



In 



A(pT < i < 1 
t < A(fT < 1 



(3.19) 



Even at vanishing temperatures, a flux dependent contri- 
bution to the persistent current is oc In(l/Aiy9) remains. 
Close to the critical mean field line (see Fig. [T]) there 
is a parametrically large enhancement of the persistent 
current due to quantum fluctuations that decays slowly 
when moving away from that line. 

In order to put this result into perspective, it is in- 
structive to compare to the persistent current in normal 
metal rings . The magnitude of the normal persistent cur- 
rent i^SEH 



IN 



1 D 1 
(j)Q B? log 5' 



(3.20) 



The asymptotic behavior of h for t <^ Atpx ^ 1 implies, 
that the persistent current due to pair fluctuations near 
ipco is parametrically larger and at low T <^T^ given by 



T° I i . (I 
(po (fico R \Aip 



(3.21) 
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FIG. 3: The modulus of the persistent current i = //[T^/^o] 
(see Eqs. |3.7| and |3.2| in the vicinity of the quantum critical 
point is shown in red as a function of the reduced tempera- 
ture t — T/T^ for four different fluxes (p close to ipco ~ 0.4. 
The inset shows the mean field phase diagram indicating the 
region relevant for this plot. The persistent current decreases 
when moving away from the quantum critical point. It dis- 
plays a pronounced maximum at low but finite temperatures. 
Thermal fluctuations grow with increasing temperature, but 
for fixed fiux the system moves further away from the criti- 
cal line. At vanishing temperatures the persistent current is 
still large and entirely caused by quantum fluctuations. Also 
shown is is (of Eq. ( |3.8[ )) in blue, which provides a very good 
approximation. The small difference between i and is is ins 
(Eq. BIO I . In this temperature regime i„s can be obtained 
from Eq. (3.241. The green line was calculated from the ap- 
proximate result of Eq. (3.151, which provides a handy esti- 



mate for the magnitude of i and the position of the maximum 
near the quantum critical point. 



where Aip — {tp — ^co)/y^cO measures the distance to the 
critical flux ipco- Since — ^/R is a number of order 

1 and = for a weakly disordered superconduc- 
tor, we find an enhancement factor of log((7) log(l/Aiy9). 
When increasing the flux at flxed temperature the persis- 
tent current decays logarithmically away from the transi- 
tion. It should be kept in mind, however, that the persis- 
tent current vanishes at ip — 0.5 due to symmetry reasons 
(see the discussion below Eq. 
approximations leading to Eq. 



3.7[) The validity of the 
3.17| is restricted to small 
AifiT, for larger AipT the full expression for is should be 
used (see Fig. |3|. 

Turning to flnite temperatures next, it is important 
that AifT is T-dependent itself and therefore, in order 
to reveal the full T-dependence of is, one should first 
find the transition line adT). We display the persistent 
current near the quantum critical point in Fig. [3j also 
comparing the different approximations and showing the 
contribution of the classical zero frequency part in the 
sum of Eq. (3.161. 



The maximum of | i | at finite T is a result of two com- 
peting mechanisms. As T grows from zero, thermal fluc- 
tuations become stronger. At the same time the distance 



to the critical line becomes larger for fixed fixed ip, which 
eventually leads to a decrease of With the help of the 
following analytic approximation, 

^l{T)^^l{l-2jlty3) (r«TO), (3.22) 

one can obtain an estimate for the position of the maxi- 
mum ipm {T) in the ^p — T phase diagram, 

ipm{T)-ifcQ_(^ A T^T— TT^A 1 2,2 



</'cO 



(3.23) 



In appendix [B] it is shown that at low temperatures 
T the nonsingular contribution to the persistent 

current ins can be written as 

Ins = 2nt}_^ dx ^^^^^ (3.24) 

where H{x,(p) = sin(27riy9)/(cosh(27ra:) — cos(27r(/3)). One 
can perform the integration in uj at zero temperature 



dy H{y/y,ip) 



X (1 + (2/7r) arctan [\n{y / ip%) / tt] ) (3.25) 



It is obvious, that is and ins have opposite signs. For a 
comparison with the singular contribution near the crit- 
ical point it is instructive to calculate the zero temper- 
ature value of ins right at the critical fiux, insi^cQ — 
0.45) — 0.17 and insiVcO — 0.4) — 0.43. Since ins is also 
a monotonously decreasing function of (p and vanishes at 
ip = Q.b we conclude that it is numerically small for all 
fluxes of our interest, \ins \ ^ The same remains true 
at flnite temperatures, see Figs. [3] |4] 

We display the temperature dependence of the persis- 
tent current in the entire temperature interval < T < 
TO in Fig. |4] For comparison, the thermal w = con- 
tribution is also shown. One can see, that at low tem- 
peratures nonzero Matsubara frequencies give a sizable 
contribution to the persistent current. When increasing 
the temperature the thermal a; = contribution becomes 
increasingly important. In Fig.|4]we compare the numer- 
ically obtained current i to the approximation is- Obvi- 
ously, it provides a very good approximation in the entire 
temperature range. 

In summary, in this section we analyzed the persistent 
current near the flux-tuned quantum critical point and 
then extended the discussion to the temperature regime 
< T < T^. Even at vanishingly small temperatures 
quantum fluctuations lead to a persistent current that is 
parametrically larger than the normal persistent current. 
When increasing the flux starting from Lp^ the persistent 
current decreases slowly and vanishes for 1^9 = 0.5. At 
this point the contributions from all angular momentum 
modes precisely cancel. When increasing the tempera- 
ture from zero at fixed fiux the persistent current in- 
creases initially since thermal fluctuations set in. This 
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FIG. 4: The modulus of the persistent current i = //[Tc/^o] 
(calculated numerically according to Eq. 3.7 in combination 



with Eq. 3.21 is shown in red as a function of the reduced 
temperature t — T/T^ for three different fluxes ip close to 
V5cO = 0.4, LP = 0.4005, 0.401, 0.402 0.41 and 0.45 from top 
to bottom. The inset shows the mean field phase diagram 
indicating the region relevant for this plot. Also displayed 
is the thermal (oj = 0) contribution to the current, which 
goes to zero at vanishing temperatures (black dashed line). 
Thermal fluctuations become increasingly important as the 
temperature grows. For two examples we indicate the singu- 
lar contribution is (of Eq. | |3.8[ )) in blue, which gives a very 
good approximation for all parameters. The small difference 
between i and is i^s- 



can be particularly well seen from Fig. |4] (dashed line) . 
On the other hand, the distance to the critical line in- 
creases for a fixed flux when increasing the temperature, 
which reduces fluctuations. The competition of these two 
mechanisms leads to a maximum in the persistent cur- 
rent at finite temperature. This maximum is the more 
pronounced the closer the flux is to the critical flux Lpc- 
We derived an approximate expression for the position 
of the maximum in the ip — T phase diagram, Eq. (3.23 1. 



It can be expected that the general features, such as the 
saturation at low T and the appearance of a maximum at 
flnite T remain intact for non-ideal rings, e.g., if the rings 
have a flnite width, because they are mainly determined 
by a single dominant angular momentum mode. In con- 
trast, the critical flux tpc may vary (see next section) and 
the sign change of the current may occur at a flux that 
is not equal to a half-integer (since many modes are in- 
volved and the phase diagram is not perfectly periodic in 
(p any more). 



IV. DISCUSSION 



persistent current in the normal state. Ref. suggests 
that a similar mechanism due to magnetic impurities is 
related to the u nexpe ctedly large persistent current in 
noble metal rings.'2SI33 

As mentioned previously, the maximal reduction of Tc 
at finite flux in the experiment of Ref. [7] was about 6 
%. It has so far not been possible to measure the persis- 
tent current close to the quantum phase transition. For 
this type of experiment one would need both sufficiently 
small rings [in order to fulfill the condition r < 0.6] and 
a measurement device that allows measuring the persis- 
tent current at the comparatively strong magnetic fields 
necessary to generate a flux of (p ~ 4>o/'^ threading the 
small area ttR^. (In the experiment at Yal^^H the use 
of cantilever required a strong magnetic field that most 
probably reduces almost totally the contributions of pair 
fiuct nations.) 

We suggest two possible strategies to relax these condi- 
tions. If the size of the ring is the main problem, one can 
try to use wider rings, because in this case the magnetic 
field penetrating the annulus of the ring helps to suppress 
superconductivity. A first consequence is that the criti- 
cal flux (fico is reduced and correspondingly the condition 
(pco < 0.5 for observing the quantum phase transition in 
the ring geometry is less restrictive, the ring radius R 
is allowed to be larger. This effect, however, is rather 
small, as we will see below. A second effect is that the 
phase diagram is no longer periodic. It can then be ad- 
vantageous to consider the transition at (/) ss l-50o or 
even higher fluxes, because then the effect of the mag- 
netic field itself (as opposed to the flux) is stronger (see 
Fig. [5] below). This approach requires measurements at 
high magnetic fields. If in turn the main problem lies 
in measuring at high magnetic fields, then the addition 
of magnetic impurities can help. Magnetic impurities re- 
duce T° itself and thereby also reduce ipco- We will now 
discuss the two mentioned effects in more detail. 

Rings of finite width: So far we considered the idealized 
case for which the width w of the ring (in radial direction) 
is vanishingly small. Next we discuss corrections to this 
result, resulting from a finite width. While doing so, we 
will still assume that the width is much smaller than both 
the coherence length ^ and the penetration depth A. The 
first assumption implies that the order parameter field 
does not vary appreciably as a function of the radius r, 
the second assumption implies that the magnetic field is 
almost constant as a function of r. 

For a ring of finite width the pair-breaking parameter 



acquires a correction, ao,, 



a, 



D 



(1) 
On ■ 



Here, a[,°j 



2]^{n — ^) is the width-independent part used so far 
and the leading corrections in w/R ar^^^ 



Our results are obtained for the case when the flux 
acts as a pair-breaking mechanism. Other pair-breaking 
mechanisms, e.g. magnetic impurities or a magnetic fleld 
penetrating the ring itself will lead to similar results. 
They cause a reduction of to zero, the pair fluctua- 
tions, however, lead to a parametric enhancement of the 



a, 



(1) 

On 



D 



w 



2i?2 4i?2 





'1 


2 1 






3^ 







(4.1) 



The relation = ttD/{8T^) can be used to express the 
result through and r. Most importantly, the width- 
dependent correction a^^J^ is not a function of n - 
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Correspondingly, the phase diagram is no longer periodic 
in if. The interpretation of this result is simple. Since 
superconductivity is already weakened by the magnetic 
field, Tc can be suppressed at a smaller flux compared to 
a ring of vanishing width. For the sake of brevity, we will 
write only the leading correction in the following. 

Let us examine some of the consequences. For the 
transition line we should now solve an equation analogous 
to Eq. (3.4), where now aQ{Lp) should be replaced by 
Q!o,m = rnin„(Q;o„). Let us first consider the regime of 
small suppression r ^ 1 and small fluxes <^ \. Then 
we can use Tc{ip) ~ T°, ao,m = cto,o and for aofi{(p) <C 
Tc{ip) one can approximately calculate the reduction of 
the transition temperature. At vanishing flux, there is 
no magnetic field and Tc{ip — 0) = is unchanged, at 
small but finite flux, however, there is a correction. 



T,M«T° 1-^ 1 



w 

Am 



(4.2) 



For small fluxes the Tc reduction is slightly stronger than 
for vanishing width. 

The condition for a suppression of T° to zero close to 
~ 1/2 can also be found. As for the case w = the crit- 
ical value for a, for which Tc vanishes, is aco = t^T^ /2^e- 
The critical flux, however, should now be found by equat- 

ing a,o to ao,o{v) = f (l + ^). The result for 
the critical flux is 



0.83 r 1 



8R^J 



(4.3) 



As expected, for a ring of finite width the critical flux is 
reduced. 

In Fig. [5] we show the mean field transition line for 
two rings with the same radius r = 0.66, but different 
widths. The ring with vanishingly small width has a flux- 
periodic phase diagram and does not show a quantum 
phase transition, since r > 0.6. The other ring has a 
width of w = i?/3. One observes three main changes. 
The maxima in Tc at finite flux are reduced compared to 

. They are shifted towards smaller flux, i.e. they do 
no longer occur at integer values of (f. Finally, the ring 
exhibits a quantum phase transition close to (p = 1.5. 

Magnetic impurities: In this article we have so far dis- 
cussed the role of an external magnetic field as the ori- 
gin of the pair-breaking mechanisms. In principle, there 
are other effects that may cause pair-breaking. Among 
them are proximity effect, exchange field, magnetic im- 
purities or interaction with the electromagnetic environ- 
ment. Each pair-breaking mechanism will have its own 
pair-breaking parameter, and to a good approximation 
the total pair breaking parameter ao-tot is the sum of 
the individual ones. The effects of these pair-breaking 
mechanisms can be obtained formally by substituting ao 
by ao-tot in the formulas discussed above. In particu- 
lar, transition temperature Tc{ip) can be obtained from 
Eq. 3.4 after the replacement Q!o(v) ~^ Q^o, tot ('/')■ The 




(f 2.0 



FIG. 5: The mean field transition line for a ring with finite 
width w = R/3 and r = 2/3 (dashed blue line), for a ring 
with vanishingly small width w = 0, but same radius r — 2/3 
(solid red line), and for a ring with vanishing width w — and 
same radius r = 2/3 in the presence of magnetic impurities 
[l/T^Ts = 0.25] (black dotted line). It can be seen that for 
the ring with finite width a quantum phase transition occurs 
near (p = 1.5, while for the ideal one-dimensional ring there is 
no quantum critical point due to the periodicity of the phase 
diagram. In the presence of magnetic impurities the phase 
diagram for the ring with w — remains periodic, but now 
quantum transitions can be found close to ip — n + 1/2 for 
any integer n. 
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FIG. 6: The mean field transition temperature of rings 
with r — 2/3 is plotted as a function of the width w/R for 
two fiuxes, p = 1.5 (black solid lines), and p = 0.5 (red 
dashed lines), and for a varying spin scattering rate 1/T^Ts 
[1/T°T, = for (a), (d), 1/T°r, = 0.1 for (&), (e), and 
l/T°Ts = 0.15 for (c), (/)]. Compared to the case p = 0.5, 
the width dependence is much stronger for p = 1.5. The ad- 
dition of a small amount of magnetic impurities reduces the 
minimal width, for which Tc vanishes. It is noteworthy that 
for a larger scattering rate than considered in this figure, a 
quantum transition can be induced near p — 0.5 even in the 
limit of vanishing width (see Fig. [5| 



condition for the quantum critical point reads ao tot{p) — 

Of particular interest is the case of magnetic impu- 
rities , disc ussed by Bary-Soroker, Entin-Wohlman and 
jj^PylHiiIl j-Qj. .^vj^ich Q;o,mi = l/'T's is cqual to the scat- 
tering rate caused by the magnetic impurities and inde- 
pendent of the flux. A sufficiently large concentration 
of magnetic impurities will reduce Tc to zero [even at 
vanishing flux] and a large persistent current is obtained 
due to the pair fluctuations. Ref. [T3] suggest that the 
large persistent current observed in copper ringj^fil may 
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Tc 



becomes 



FIG. 7: The modulus of the susceptibUity for — 1.5 and 
T = 0.5 r^" is plotted as a function of the width w/R for three 
values of the spin scattering rate, 1/T^Ts = 0.0 for (a), 0.1 for 
(6), and 0.2 for (c). The radius of the ring is r = 2/3. Also 
shown is the width-dependence of the mean field transition 
temperature Tc/T^ at ip = 1.5 (dashed lines) for two values of 
the spin scattering rate, I/T^Ts = for (a), and = 0.1 

for For l/T^r^ = 0.2, Tc is equal to zero for arbitrary 
width. The susceptibility at ^ = 1.5 is entirely caused by 
fluctuations, because the mean field transition temperature is 
smaller than 0.5Tc , even in the limit of vanishing width. We 
see that the susceptibility is a smooth function of the width. 
In particular, it remains large even after passing the threshold 
values of w/R for which Tc vanishes. 



be attributed to such pairing fluctuations. 

As emphasize earlier, in our case the addition of mag- 
netic impurities may push the system to the quantum 
critical point at a smaller external magnetic field, since 
it provides an additional flux-independent pair-breaking 
mechanism. If the main experimental difliculty is to per- 
form sensitive measurements at high magnetic fields, in- 
troducing magnetic impurities might therefore make pos- 
sible the experimental observation of the flux tuned quan- 
tum critical point, see Fig. |6] 
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Appendix A 

In this appendix we will show how to calculate the 
partition function when two modes "00 and tpi are taken 
into consideration. We start from Eq. (2.8 1 and introduce 
the notation ^„ = a„/r°, B = b/{2VTY), x = |V'oP and 



Z = tt" 



dxdy e 



-Aox-Aiy-B{x'^+y^+'kxy) 



(Al) 



We want to perform one integration explicitly. To this 
end we first change integration variables to w = x + y and 
z = x — y and obtain Z — dw J^^^ dz exp(— wA^ — 

zA_ - -f iSz^), where A± = {Aq ± yli)/2. After 

changing the order of integration we find 



Z = 



dz 



(A2) 



xe"^' cosh(^_zV2AB) erfc(V32 + A+/^/QB) 



Next we derive a formula for persistent current. Combin- 
ing the formulas / = (2mi?2(^o)"^ Y.n=o,i{'^- v) 
and (IV'nP) = ^^Ar, InZ one finds 



1 



2A(^9. )lnZ, 



AmR^Tl 

where = (/? - 1/2 = AmR'^T^A^. Let us note that 
A 



(A3) 



dA+ In Z = 
dA_ InZ = 



SB 

A_ 



UB^/^Z 



e''oerfc(a:o) + e^ierfc(xi) 

e'^oerfc(a;o) — e^ierfc(xi) 
(A4) 



where x„ — An/{2\/B) = Sn/Gi coincides with the vari- 
able Xn used in the main text. After combining these 
results one finds 



12 = -4— Al_ + -X-M- 
Gi 3 



where 
M± 

and 



± -r .p 



1/2 2 

- l^e^^oeric^xo) ± e^i erfc(a;i) 



(A5) 



(A6) 



POO 

Jo 



z + xi) (A7) 



where x± = (a;oia^i)/2. An analogous formula has been 
given in Ref . HI 

Most interesting for us is the quantity x(l/2) = 
— di2/d(p\^^-j^^2- It is worth noting that for = 1/2 
one finds A- = x- = (i.e xq = Xi), and formulas 
simplify considerably. We can obtain the formula for 
x(l/2) directly by differentiating the result for Z2 or by 
first differentiating Eq. ( A3 ) and then using the relations 



in Eq. (A4). In the latter case one obtains as an in- 

J 2 

termediate step the relation — — 2mR°T°z (^^+ 

8T°mR^ ) ^\a-=o- result for x is stated in the 



y = iV'il ■ Then the expression for the partition function main text, Eq. (2.17) 
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Appendix B 

In this appendix we sketch the derivation of the ex- 
pressions for the persistent current in the Gaussian ap- 
proximation used in the main text. We also make con- 
tact with Ref. 1151 where fluctuations at temperatures 
T > Tc{(p = 0) have been examined [In the presence of 
magnetic impurities, Tc{ip = 0) differs from T^?, the tran- 
sition temperature in the absence of any pair-breaking 
mechanism, and may even vanish] . Our starting point is 
Eq. ISJl), 



^ £-i(n,w) 



(Bl) 



Using standard methods for transforming the sum in n 
into an integral in the complex plane, we can write i as 



iG = — dz H{z,tp) R{z,u}) 



(B2) 



where 



R{z,uj) 

H{z,^) 
and we use the notation 



dzC~^{—iz, uj) 
C~^{—iz,iij) 
sin(27r(p) 
cosh(27rz) — cos(27riy9) 



(lyC) ^{—iz,u}) = 



(B3) 



(B4) 



In 



■ Etz 



AttT 



The contour C' includes all the poles of H{z, ip) with 
'<sz > 0. These poles are located at 



z = i(p,i{l - <^),i(l + ip),i{2-tp),.. 



(B5) 



The contour C' does not include poles of R{z,uj), see 
Fig. [8j Our goal will be to deform the contour of inte- 
gration in such a way that the integral can be evaluated 
at the poles of R with the help of the residue theorem. 
Poles of R{z,uj) occur either due to zeros or poles of 
C^^{—iz,Ld). Due to the dependence of on the 
poles of R{z, w) come in pairs, for each pole at a point z 
there is a pole at —z. 

Zeros of £^^{—iz, cj) can be classified according to the 
value of 



A{z,u) ^ 1/2 + {\uj\ - eTZ^)/{47rT), 



(B6) 



the argument of the digamma function. There is a pair 
of zeros ±Zs, for which A is positive, and there is a pair 
of zeros ±z„ for each interval —n ~ 1 < A < —n, where 
n = 0,1,2,.... The existence of such zeros is clear from 
the reflection formula 



-0(1 — z) — ip(z) + 7rcot(7rz). 



(B7) 



i{l-Hp) 



c 



FIG. 8: The complex z plane. The poles ol H{z,ip) (Eq. B3| 



as a function of z are indicated by crosses together with the 
orig inal integration contour C' . Equivalently, the integral in 
Eq. B2 can be calculated by integrating along the path C" , 
and additionally taking into account the residue at the pole 
Zs- Zs can be either imaginary or real and is responsible for 
the dominant contribution is to the persistent current. The 
contour C" encloses Zn and S„, the poles of R{z,uj), which 
are positioned on the real axis to the right of zq. As the 
temperature decreases the distance between adjacent z^ and 
Zn becomes small and the sum over the corresponding residues 
can be approximated by an integral along a branch cut. 



Poles of £ ^ {^iz, uj) originate from poles of the digamma 
function for A = — n, n = 0, 1, 2, . . . and we label the pairs 



accordingly, ±z„, where 



= ^({2n + l)2TiT+\Lo\)leT. 



(B8) 



Note that all z„ and 2„ are real and z„ < z„ < 

The zeros at iz^ deserve special attention. For T > 
r°, Zs is real. For T < T^, however, Zg can be ei- 
ther real or imaginary. Indeed, due to the relation 
ln(r/T°) = V-(l/2)-^(l/2 + ac{T)/2TTT), we can write 
Zs = (pc{T)y/\uj\/2a{T) ~ 1 for T < T° . Zs is real for uj > 
2a{T) and purely imaginary Zg = i^pc{T)^yl — \uj\/2a{T) 
for UJ < 2a{T). 

By deforming the contour C and using the symmetry 
and convergence-properties of the integrand we can write 
is as an integral along a contour C" that encloses all poles 
of R{z, uj) on the positive imaginary axis for < '<sz < ip, 
and on the positive real axis in the complex z-plane, see 
Fig. [8) The integral can then formally be calculated with 
the help of the residue theorem. 

When UJ is small and Zs imaginary, H can become large 
for (fi ^ (fis, reflecting the closeness to the phase transi- 
tion. Evaluating the residue at Zs one finds the singular 
contribution to the persistent current 

Is = -27rtJ2H{zs,ip) (B9) 
The nonsingular contribution arises from all other poles 



15 



of R. and can be written as 



^-2TTtJ2Y.^H{z„,^)-HiS„,^)]. (BIO) 



In general, the poles z„ need to be determined nu- 
merically. The derived representation is particularly use- 
ful in two limiting cases, either for very high or very 
low temperatures. For very high temperatures Zs is real 
the spacing between consecutive z^, z„, z„ is large. Since 
H{z, if) decays fast as a function of z one can approxi- 
mate the result well by considering just the smallest of 
the and in this way obtain relatively simple 

formulas. This approximation was utilized in Ref. 1151 
where a similar representation was derived for tempera- 
tures T > Tc{ip = 0) (in the presence of magnetic impu- 
rities), when all poles of R{z,uj) are real. Another useful 
limit is the limit of very low temperatures, when Zs is pos- 
sibly imaginary, but the other z„, z„ are closely placed 
on the real axis. Then, calculating the residues for i„s is 
very similar to an integration along a branch cut as we 
will describe now. In this low temperature limit one can 
write 



Simple substitution gives the formula Eq. (3.24) stated 
in the main text. In the limit T — > 0, it is more 
convenient to use the relation [y([n^ y + tt^)]^^ = 
arctan (ln(j/)/7r), perform a partial integration in y, 
subsequently use dyf{y+\uj\/2ac{T)) = 2ac{T)di^if{y + 
\uj\/2ac(T)) to perform the integral in a; and combine the 
result with the boundary term. The result is Eq. (3.25). 

Finally, let us make contact with the analysis of 
Ref. [13 In this paper, fluctuations were analyzed for high 
temperatures T > Td f = 0) (in the presence of magnetic 
impurities) and the flux harmonics of the persistent 
current i = im sin(27rrniy9) were calculated. In this 
situation one may use the Poisson summation formula to 
perform the sum over angular momentum modes, since 
the flux dependence of i is smooth for T ^ Tc(</? = 0). 
[This is not so for T < Tc{(p = 0) due to the phase transi- 
tion to the superconducting state, or, in a mathematical 
language, due to the presence of the pole at purely imag- 
inary Zs-] A comparable situation arises for T > in 
the absence of magnetic impurities: All poles of R{z,uj) 
(including Zg) are real, one can use Eqs. B9 and BIO and 
the relation 



« In 



1 

2 ' 
an 



an + 



2ttT 



H 

2 



ac{T) 



2ttT 



(Bll) 



H{x, ip) = 2j2 e^^'r:!™ sin(27rm(/7) (B13) 



which is valid for a; > 0, to write 



Following the same analysis presented above with this 
approximate representation, is remains unchanged. For 
ins we have to perform an integration along a branch cut 
appearing for z > zq ^ ^/\uJ\JeT, which is related to the 
branch cut of the logarithm. It gives the result 



dy H{ip,{T)^y+\u\/2a,{T ),ip\ 
y In^ 2/ + 



i = —Ant sin(27rm(p) 



(B14) 



E 



{B12) from which one can read off the flux harmonics 
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